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Abstract 

The axino and the gravitino are weh-motivated candidates for the hghtest supersymmetric par- 
ticle (LSP) and also for cold dark matter in the Universe. Assuming that a charged slepton is the 
next-to-lightest supersymmetric particle (NLSP), we show how the NLSP decays can be used to 
probe the axino LSP scenario in hadronic axion models as well as the gravitino LSP scenario at 
the Large Hadron Collider and the International Linear Collider. We show how one can identify 
experimentally the scenario realized in nature. In the case of the axino LSP, the NLSP decays will 
allow one to estimate the value of the axino mass and the Peccei-Quinn scale. 
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I. INTRODUCTION 



In supersymmetric extensions of the Standard Model with unbroken R-parity [1], the 
hghtest supersymmetric particle (LSP) is stable and plays an important role in both collider 
phenomenology and cosmology. The most popular LSP candidate is the lightest neutralino, 
which appears already in the Minimal Supersymmetric Standard Model (MSSM). Here we 
consider two well-motivated alternative LSP candidates, which are not part of the spectrum 
of the MSSM: the axino and the gravitino. In particular, either of them could provide the 
right amount of cold dark matter in the Universe if heavier than about 1 MeV (see [2, 3] 
and [4-8], respectively, and references therein). 

The axino [9-11] appears (as the spin-1/2 superpartner of the axion) when extending 
the MSSM with the Peccei-Quinn mechanism [12] in order to solve the strong CP problem. 
Depending on the model and the supersymmetry (SUSY) breaking scheme, the mass of the 
axino can range between the eV and the GeV scale [10, 13-15]. 

The gravitino appears (as the spin-3/2 superpartner of the graviton) once SUSY is pro- 
moted from a global to a local symmetry leading to supergravity (SUGRA) [16]. The mass 
of the gravitino depends strongly on the SUSY-breaking scheme and can range from the 
eV scale to scales beyond the TeV region [1, 17, 18]. In particular, in gauge- mediated 
SUSY breaking schemes [17], the gravitino mass is typically less than 100 MeV, while in 
gravity-mediated schemes [1] it is expected to be in the GeV to TeV range. 

Both the axino and the gravitino are singlets with respect to the gauge groups of the 
Standard Model. Both interact extremely weakly as their interactions are suppressed by 
the Peccei-Quinn scale [19, 20] /„ > 5 x lO'' GeV and the (reduced) Planck scale [20] Mpi = 
2.4 X 10^^ GeV, respectively. Therefore, in both the axino LSP and the gravitino LSP 
cases, the next-to-lightest supersymmetric particle (NLSP) typically has a long lifetime. For 
example, for axino cold dark matter, an NLSP with a mass of 100 GeV has a lifetime of 
0{1 sec). For gravitino cold dark matter, this lifetime is of 0{1 sec) for a gravitino mass 
of 10 MeV and of O{10^ sec) for a gravitino mass of 10 GeV. Late NLSP decays can spoil 
successful predictions of primordial nucleosynthesis and can distort the CMB blackbody 
spectrum. Constraints are obtained in order to avoid the corresponding (rather mild) axino 
problem or the more severe and better-known gravitino problem. In the axino LSP case, 
either a neutralino or a slepton could be the NLSP [21]. In the gravitino LSP case, these 
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constraints strongly disfavour a bino-dominated neutralino NLSP, while a slepton NLSP 
remains allowed [8, 22]. 

Because of their extremely weak interactions, the direct detection of axinos and gravitinos 
seems hopeless. Likewise, their direct production at colliders is very strongly suppressed. 
Instead, one expects a large sample of NLSPs from pair production or cascade decays of 
heavier superparticles, provided the NLSP belongs to the MSSM spectrum. These NLSPs 
will appear as quasi-stable particles, which will eventually decay into the axino/gravitino 
LSP. A significant fraction of these NLSP decays will take place outside the detector and will 
thus escape detection. For the charged slepton NLSP scenario, however, there have recently 
been proposals, which discuss the way such NLSPs could be stopped and collected for an 
analysis of their decays into the LSP. It was found that up to 0(10^-10'^) and (9(10^-10^) 
of charged NLSPs can be trapped per year at the Large Hadron Collider (LHC) and the 
International Linear CoUider (ILC), respectively, by placing 1-10 kt of massive additional 
material around planned collider detectors [23, 24]. 

In this Letter we assume that the NLSP is a charged slepton. In Sec. II we investigate the 
NLSP decays in the axino LSP scenario. These decays were previously considered in [21]. 
We show that the NLSP decays can be used to estimate the axino mass and to probe the 
Peccei-Quinn sector. In particular, we obtain a new method to measure the Peccei-Quinn 
scale fa at future coUiders. 

In Sec. Ill we consider the corresponding NLSP decays in the gravitino LSP scenario. 
These decays were already studied in [25] . It was shown that the measurement of the NLSP 
lifetime can probe the gravitino mass and can lead to a new (microscopic) determination 
of the Planck scale with an independent kinematical reconstruction of the gravitino mass. 
Moreover, it was demonstrated that slepton NLSP decays into the corresponding lepton, the 
gravitino, and the photon can be used to reveal the peculiar couplings and possibly even the 
spin of the gravitino. In Ref. [25] the limit of an infinite neutralino mass was used. Here we 
generalize the result obtained therein for the three-body decay by taking into account finite 
values of the neutrahno mass. 

A question arises as to whether one can distinguish between the axino LSP and the 
gravitino LSP scenarios at colliders. From the NLSP lifetime alone, such a distinction will 
be difficult, in particular if the mass of the LSP cannot be determined. Thus, an analysis of 
the three-body decay of the charged NLSP slepton into the corresponding lepton, the LSP, 
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and a photon will be essential. With a measurement of the polarizations of the final-state 
lepton and photon, the determination of the spin of the LSP should be possible [25] and 
would allow us to decide clearly between the spin- 1/2 axino and the spin-3/2 gravitino. 
The spin measurement, however, will be very difficult. In Sec. IV we present more feasible 
methods to distinguish between the axino LSP and the gravitino LSP scenarios, which are 
also based on the analysis of the three-body NLSP decay with a lepton and a photon in the 
final state. 

Let us comment on the mass hierarchy of the relevant particles. There are six possible 
orderings in the hierarchy of the axino mass m~, the gravitino mass mg, and the mass 
of the lightest ordinary supersymmctric particle (LOSP) mLosp- Here the LOSP is the 
lightest charged slepton. The cases relevant in this Letter are (i) m~ < mLosp < ''^g; 
(n) mg < mLosp < rn-, (in) m~ < uiq < mLosp, and (iv) mg < m~ < mLosp- In cases 
(iii) and (iv), the LOSP has two distinct decay channels, one into the axino and the other 
into the gravitino. However, unless the decay rates into the axino and the gravitino are 
(accidentally) comparable, the phenomenology of the LOSP decay in the cases (iii) and (iv) 
can essentially be reduced to the cases (i) or (ii), although not necessarily respectively, as 
will be discussed in Sec. IV. We will thus concentrate on the cases (i) and (ii) and call the 
LOSP the NLSP. 

II. AXINO LSP SCENARIO 

In this section we consider the axino LSP scenario. The relevant interactions of the axino 
are discussed. The rates of the two-body and three-body decays of the charged slepton NLSP 
are given. We demonstrate that these decays can be used to estimate the Peccei-Quinn scale 
and the axino mass. 

To be specific, we focus on the case where the lighter stau r is the NLSP. In general, 
the stau is a linear combination of tr and tl, which are the superpartners of the right- 
handed and left-handed tau lepton, respectively: f = cos^t-Tr -I- sin Ot-Ti,. For simplicity, 
we concentrate on a pure 'right-handed' stau tr, which is a good approximation at least 
for small tan j3. Then, the neutralino-stau coupling is dominated by the bino coupling. In 
addition, we assume for simplicity that the lightest neutralino is a pure bino. 
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A. Axino Interactions 



Let us first discuss how the axino couples to the stau. Concentrating on hadronic, or 
KSVZ, axion models [26] in a SUSY setting, the coupling of the axino to the bino and 
the photon/Z-boson at scales below the Peccei-Quinn scale is given effectively by the 
Lagrangian [2] 

>C~ = i ^^Q^J^ « 75 bfM: lu] B (cos OwFt^i, - sin 6wZ,j,i,) , (1) 

where 9w is the weak mixing angle, aY — o./ cos^ 9w with the fine structure constant a, 
and F^j, and Z^y are the field strength tensors of the photon and Z-boson, respectively. The 
interaction Lagrangian (1) is obtained by integrating out the heavy (s) quarks introduced 
in super symmetric KSVZ axion models. Indeed, the KSVZ axino couples directly only 
to these additional heavy (s) quarks. Thus, the above coupling depends, for example, on 
the hypercharge of these heavy (s)quarks, which we assume to be non-zero. The model 
dependence related to the Peccei-Quinn sector is expressed in terms of the factor CaYY — 
As the MSSM fields do not carry Peccei-Quinn charges, the axino couples to the stau 
only indirectly, via the exchange of intermediate gauge bosons and gauginos. 

In the alternative DFSZ axion models [27], once supersymmetrized, the mixing of the 
axino with the MSSM neutralinos can be non-negligible and other couplings between the 
axino and the MSSM fields will arise. Here, however, we focus on the KSVZ-type models. 



B. The Two-Body Decay t ^ t -\- a 

We now consider the two-body decay f ^ t -\- a m. the framework described above. 
We neglect the tau mass for simplicity. With the effective vertex (1), i.e. with the heavy 
KSVZ (s)quarks integrated out, this two-body decay occurs at the one-loop level. The 
corresponding Feynman diagrams are shown in Fig. 1, where the effective vertex is indicated 
by a thick dot. Using the method described in [28], we obtain the following estimate for the 
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FIG. 1: The dominant contributions to the two-body NLSP decay tr — > r + a. 
decay rate:^ 

where is the mass of the bino and m~ is the mass of the stau NLSP, i.e. m~ < m~ < m^. 
As explained below, there is an uncertainty associated with the method used to derive the 
decay rate (2). We absorb this uncertainty into the mass scale m ~ — (^(100 GeV) 

and into the factor ^ ^ in the first line. We used log(/a/m) ~ 20.7 to get from the 
first to the second line. 

Here a technical comment on the loop integral is in order. If one naively integrates 
over the internal momentum in the diagrams with the effective vertex — see Fig. 1 — one 
encounters logarithmic divergencies. This is because the effective vertex (1) is applicable 
only if the momentum is smaller than the heavy (s) quark masses, whereas the momentum 
in the loop goes beyond that scale. In a rigorous treatment, one has to specify the origin 
of the effective vertex, i.e. the Peccei-Quinn sector, and to calculate the two-loop integrals 
with heavy (s)quarks in the additional loop. Such a two-loop computation leads to a finite 
result [29]. Here, instead, we have regulated the logarithmic divergencies with the cut-off 

1 We correct the factor of (1/16)(1 + tan^ 6*^0^/(1 - tan^ Owf-, which is missing in Eq. (3.12) of Ref. [21]. 
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FIG. 2: The dominant contributions to the three-body NLSP decay tr — > r + 7 + a. 

fa and kept only the dominant contribution. The mass scale m and the factor ^ have been 
introduced above to account for the uncertainty coming from this cut-off procedure. 

C. The Three-Body Decay r — ^ t + 'y + a 

We now turn to the three-body decay tr r + 7 + a. We again neglect the tau mass for 
simplicity. In contrast to the two-body decay considered above, the three-body decay occurs 
already at tree level, once the effective vertex given in (1) is used. In addition, we take into 
account photon radiation from the loop diagrams of Fig. 1, since the additional factor of 
a is partially compensated by the additional factor of log(/a/m). As above, we keep only 
the dominant contribution of the loop diagrams. The corresponding Feynman diagrams are 
shown in Fig. 2, where a thick dot represents the effective vertex (1) and a shaded triangle 
the set of triangle diagrams given in Fig. 1. As the photon radiation from an electrically 
charged particle within the loops leads to a subdominant contribution, these processes are 
not shown in Fig. 2. At each order in log(/a/m), only the leading order in a is computed 
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while higher-order corrections are not considered. In terms of the observables that seem 
to be most accessible, i.e. the photon energy and cos^, the cosine of the opening angle 
between the photon and the tau direction, the corresponding differential decay rate reads 

spins 



dx^dcosO 5127r3 [1 - (a;^/2)(l - cos 



Y.\M{f^^T^~a)\\ (4) 



where 



E |M(rR-r7S)r = -—^ ^ cos^, , (5) 

TT COb J a 



spins 



with 



X = ^ A~=^ A~='^ (6) 
m~ mA- ^ mi- 

~ T T 



and 



(jT) _ x^l-A~-x^)[l+cose+A~{l-cos9)][l+cos9+A^{l-cos9)] 

^^^[X^, COSU, A^J - ^^^^^ ^ ^^^^^ ^ _ ^^^2 _ ^^(1 _ cos^)]}2 



+ 77;- Clog 



TTCOS^^i^ \m J [x-y(l + cos^) + 2^4- — ^4^(2 — a;-y(l — cos^)] 

^ .4^[(l + cos/?)(l-A^) + A^,r,(l-cos^?)] 



+ cos 0) + 2A~ - A^[2 - x-y{l - cos 9)] 



9a' 2 .2 (IAaJ 1 + cos9 + A~{1-cos9) 2(1 + cosg)(l - 1 
47r2cos4^w'^ ^ [mj ^\ (1 - cos^)(l - ^~ - x^) ^^(l-cos^) J" 

Hereafter, wc use log (fa/m) ~ 20.7, as in the previous section. 

The three-body decay r — > r+7 + a involves bremsstrahlung processes (see Fig. 2) and, as 
already mentioned, we have neglected the tau mass. Thus, when the photon energy and/or 
the angle between the photon and the tau direction tend to zero, there are soft and/or 
coUinear divergences. Consequently, the total rate of the decay f— >T-|-7-|-ais not defined. 
We define the integrated rate of the three-body decay f — > r + j + a with a cut on the scaled 
photon energy, x^ > a;™*, and a cut on the cosine of the opening angle, cos^ < 1 — x^"*: 

r TR^T7a;4"*,a;^"* = / dx^ dcos9 ^ . 8 

' A^"* i-i dXjdcos9 

As explained in Sec. IV, the quantity r(Ta — > r ^ a ; x^^ , Xq^^) will be important in distin- 
guishing between the axino LSP and the gravitino LSP scenarios. 
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D. Probing the Peccei— Quinn Scale and the Axino Mass 



In the axino LSP scenario, the stau NLSP decays provide us with a new method to probe 
the Peccei-Quinn scale /„ at coUiders. As we will see in Sec. IV B, the branching ratio of the 
three-body decay is small if reasonable cuts are used. Thus, we can use the two-body decay 
rate (3) to estimate the stau lifetime, r- ~ l/r(T ra). Accordingly, the Peccei-Quinn 
scale fa can be estimated as 




once m~, m^, and the lifetime of the stau r~ have been measured. The dependence on the 
axino mass is negligible for m~/m~ < 0.1, so that /„ can be determined without knowing 
m~. For larger values of m~, the stau NLSP decays can be used to determine the mass of 
the axino kinematically. In the two-body decay f — > r -|- a, the axino mass can be inferred 
from the energy of the emitted tau lepton: 

^m^ + ml- 2m~£^T- , (10) 

with an error depending on the experimental uncertainty on m~ and E^-. 

III. GRAVITINO LSP SCENARIO 

In this section we assume that the gravitino is the LSP and again that the pure right- 
handed stau is the NLSP. The corresponding rates of the two-body and three-body decay 
of the stau NLSP are given. These decays have already been studied in Refs. [25] . Here we 
generalize the result obtained for the three-body decay by taking into account finite values 
of the neutralino mass. For simplicity, we assume again that the lightest neutralino is a pure 
bino. 

The couplings of the gravitino G to the tr, r, B, and 7 are given by the SUGRA 
Lagrangian [16]. The interactions of the gravitino are determined uniquely by local SUSY 
and the Planck scale and, in constrast to the axino case, are not model-dependent. 

A. The Two-Body Decay f ^ r -|- G 

In the gravitino LSP scenario, the main decay mode of the stau NLSP is the two-body 
decay f ^ r -\- G. As there is a direct stau-tau-gravitino coupling, this process occurs at 
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tree level. Neglecting the r-lepton mass rrir, one obtains the decay rate: 

r(TR ^ tG) = - — ^ 1 - ^ (11) 

487r m ~ Mpj \ m~ / 



G 



2 / ^2 \ 4 



= (5.89 sec ^ — 1 f . (12 

^ ' VlOOGcVy V / \ ^/ 

In order to get from the first to the second line, we have used the value of the reduced Planck 
mass Mpi = (SttGn)"^^^ = 2.435 x 10^^ GeV as obtained from macroscopic measurements 
of Newton's constant [20] Gn = 6.709 x 10~^^GeV~^. Thus, the gravitino mass can be 
determined once the stau NLSP lifetime governed by (12) and m~ are measured. As pointed 
out in Refs. [25], expression (11) can also be used the other way around, i.e. for a microscopic 
determination of the Planck scale once the masses of the gravitino and the stau are measured 
kinematically. Note the strong dependence on and m~. In the axino LSP scenario, the 
corresponding rate (2) becomes independent of the axino mass for m~/m~ < 0.1, so that the 
Peccei-Quinn scale can be determined even if m~ is too small to be inferred kinematically. 

B. The Three-Body Decay r ^ r + 7 + G 

Let us now turn to the three-body decay tr — > r + 7 + G. The corresponding Feynman 
diagrams are shown in Fig. 3. We neglect again the tau mass for simplicity. For finite bino 
mass, we obtain the following differential decay rate 

(fV {tyi T ^ G) ruf Xy{l — AQ — Xy) 

cos6 5127r2 [1 — (a;-j,/2)(l — cos6 ^^^^ 



^\Mifn^r^G)\', (13) 



where 



E \Mif^^rjG)\' = ^-^ Fg\x,,cose,A^,As) (14) 

spins PI G 



with the definitions of and A^ given in (6), A^ = m~/mq:, and 

F^^^(x coze A A)- ZA' IxA 2(2-5cosg)Ag a:,(l + cosg) 
t^^{x^,co^e,A^,A^)- 6A^ lx,A^+ {I -cose) 

(l + cos^)(3 + cos^) 2(1- Ag)3(l + cos^) A^(l-Ag)2 



(1 — cos^)2 a;^(l — cos^) l — A^ — x^ 

(1 - Ag)2(l + cos^) 4 [l + cos^ + Ag(l -cos^)]' 



(1 - - x^){l - cose) [2 - x^{l - cose)Y (1 - cos^)2 

2 {3 + cos ^ [4 - cos ^ + 2^g(l - cos e)] } [1 + cos ^ + Ag(l - cos e) 



[2-a;^(l-cos^)](l-cos^)^ 
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FIG. 3: The three-body NLSP decay tr ^ r + 7 + G. 



+2(1 -A; 



l + x^-x'i- 2Ag(l + - 2x'i) + A^(l + hx^) 



x^(l - A^)(l - A 



G 



1 + X^(2 + A^) - + 2ylg(l - X, 



+ 



4(1- Ag- a; 



[2 - ^^,(1 - cos e)] ' [2 - x^(l - cos e)f 

2(1 + cose)(l - A^) + 3x^Ag(l - cos^) 



x^(l + cos 6*) + 2{A^ - A^) + A^x^{l - cos 61) 
2 {A| [-3 - + A^(2 + x^)] + ^A^A^{1 + x^ 



x: 



x^(l — Aj^) 
2 

+ 


X-y{ 


1 + cos 6*) + 2{A^ - A^) + A^x-y 
(1 - x^)(l + 2Ag + x^) + x^A^ 


(1 -cos^)" 


x^(l — Aj^) 


x^(l + COS 6*) + 2{A^ - A^) + A^x^il - cos 6*) 



4(1 - - x^ 



1 + 3A~)(1-A 



G> 



(1-A 



+ 



2 — x^ — 2(yl^ — Aj^] 



B) 



x^(l — cos6') 



2{A^-A^) 


3Ag(2 


2 j/1 Q CC^y 


) + A^{2-2A^ + x,) 






x^(l + cos t 


+ 2(Ag - 


Aj^) + A^x^(l — cos 9) 



4(1 - - x^)(3Ag + A^){Ag - A^ 



(1 - A^) x^(l + cos6') + 2(v4g - v4^) + A^x^(l - cos( 



(15) 
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In the limit — > oo, only the terms in the first four lines of (15) remain and the result 
given in the appendix of the first reference in [25] is obtained. For finite values of the bino 
mass, the diagram with the bino propagator in Fig. 3 has to be taken into account, which 
then leads to our more general result. 

As in the axino case, the total rate of the three-body decay r — r + 7 + Gis not 
defined. We thus introduce again the integrated rate with a cut on the scaled photon 
energy, > x™*, and a cut on the cosine of the opening angle, cos^ < 1 — x^"*, 

r(fR^r7G;x-\xr) = f'^" dx, /''^'^" dcos /'^f^ ^ ^] ■ (16) 

' ia;^"t J-1 dx^d COS 6 

This quantity will be used in our comparison of collider signatures of the axino LSP and the 
gravitino LSP scenarios. 

IV. AXINO VS. GRAVITINO 

In this section we show how the two-body and three-body decays of the stau NLSP can 
be used to distinguish between the axino LSP scenario and the gravitino LSP scenario at 
colliders. We compare the total decay rates of the stau NLSP, the branching ratios of the 
three-body decays f — > t + ^ + a/G with cuts on the observables, and the differential 
distributions of the decay products in the three-body decays. 

A. Total Decay Rates 

Let us discuss the lifetime of the stau NLSP in the axino LSP and in the gravitino LSP 
scenarios, and examine whether the lifetime can be used to distinguish between the two. In 
both cases, the total decay rate of the stau NLSP is dominated by the two-body decay, 

rf^tix-^irn^Tt) , i = d,G, (17) 

with the rates given respectively in (3) and (12). Thus, the order of magnitude of the stau 
NLSP lifetime is (essentially) determined by m~, mg, and in the axino LSP scenario 
and by m~ and mg in the gravitino LSP scenario. Among those parameters, one should be 
able to measure the stau mass m~ and the bino mass by analysing the other processes 
occurring in the planned collider detectors. Indeed, we expect that these masses will already 
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be known when the stau NLSP decays are analysed. To be specific, we set these masses 
to m~ — 100 GeV and — 110 GeV, keeping in mind the NLSP hfetime dependencies 
T~ oc l/(m~m~) for the axino LSP and t~ oc l/m~ for the gravitino LSP. Then, the order of 
magnitude of the stau NLSP hfetime is governed by the Peccei-Quinn scale /„ in the axino 
LSP scenario and by the gravitino mass in the gravitino LSP scenario. 

In the axino LSP scenario, the stau lifetime varies from C(0.01 sec) to 0(10 h) if we 
change the Peccei-Quinn scale /„ from 5 x 10^ GeV to 5 x 10^^ GeV, as can be seen from (3). 
For the given values of m~ and mg, these values can probably be considered as the lower 
and upper bounds on the stau NLSP lifetime in the axino LSP case. 

In the gravitino LSP case, the stau lifetime can vary over a much wider range, e.g. from 
6 X 10^*^ sec to 15 years by changing the gravitino mass from 1 keV to 50 GeV, as can be 
seen from (12). Therefore, both a very short stau NLSP lifetime, r- ;$ msec, and a very long 
one, T~ > days, will point to the gravitino LSP scenario. For example, in gravity-mediated 
SUSY breaking models, the gravitino mass is typically (10-100) GeV. Then, the lifetime of 
the NLSP becomes of C (years) and points clearly to the gravitino LSP scenario. 

On the other hand, if the observed lifetime of the stau NLSP is within the range 
0(0.01 sec)-(9(10 h), it will be very difficult to distinguish between the axino LSP and 
the gravitino LSP scenarios from the stau NLSP lifetime alone. In this case, the analysis of 
the three-body NLSP decays will be crucial to distinguish between the two scenarios. 



B. Branching Ratio of the Three-Body Decay Modes 

We now consider the branching ratio of the integrated rate of the three-body decay 
f^T + 'j + a/G with cuts 



BRif^ ^r^^■x-\ x^) ^ "^^"^^ ^ IZ ' ^^^^ ^ ^ = S, G . (18) 



In Fig. 4 this quantity is shown for the gravitino LSP (sohd fine) and the axino LSP (dashed 
line) for m~ = 100 GeV, mg = 110 GeV, = 10" GeV, ^^C^yy = 1. m^/m^ < 1, and 
mg = lOMeV.^ In the left (right) part of the figure we fix x"^^ = 0.1 (a;^"* = 0.1) and vary 



^ The results shown in Fig. 4 are basically independent of the Peccei-Quinn scale /„ and the gravitino mass 
m~ provided m~/m~ < 0.1. For larger values of the gravitino mass, the stau NLSP lifetime being of 
O(years) points already to the gravitino LSP scenario as discussed above. 
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Branching Ratios of tr t j a/G with Cuts 




FIG. 4: The branching ratio of the integrated rate of the three-body decay r^r + 7 + a/G with 
cuts as a function of x^"* for x^"* =0.1 (left) and as a function of for x^"* = 0.1 (right). The 
soUd and dashed hues show the results for the gravitino LSP and the axino LSP, respectively, as 
obtained with m~ = 100 GeV, m~ = 110 GeV, /„ = 10^^ GeV, C^Cfyy = 1, < 1, and 

mp; = lOMeV. 

Xg"* (x™*). The dependence of the branching ratio (18) on the cut parameters in the axino 
LSP case differs quaUtatively from the one in the gravitino LSP case. Moreover, there is 
a significant excess of BR{tr —>■ rja; x"^^, x™*) over BR{t}^ tjG ; x™*, x^"*) over large 
ranges in the cut parameters. For example, if 10^ stau NLSP decays can be analysed and 
the cuts are set to = 0.1, we expect about 165±13 (stat.) tr ^ T7a events for 

the axino LSP and about lOOilO (stat.) tr — > rjG events for the gravitino LSP. Thus, 
the measurement of the branching ratio (18) would allow a distinction to be made between 
the axino LSP and the gravitino LSP scenarios. For a smaller number of analysed stau 
NLSP decays, this distinction becomes more difficult. In addition to the statistical errors, 
details of the detectors and of the additional massive material needed to stop the staus and 
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Axino LSP Scenario Gravitino LSP Scenario 

1 d^rCxa^T-yg) 1 d^r(TR->T -y G) 

r(fR— ►TTo; a;="*=:c§"*=0.1) dx^dcosO r(-fH^T -y G; a;="*=xg"*=0.1) dx^dcos9 




FIG. 5: The normalized differential distributions of the visible decay products in the decays r r+ 
^ + a/G for the axino LSP scenario (left) and the gravitino LSP scenario (right) for m~ = 100 GeV, 
= 110 GeV, rn~/rn~ -C 1, and m~ = lOMeV. The cut parameters are set to x™* = x™* = 0.1. 
The contour lines represent the values 0.2, 0.4, 0.6, 0.8, and 1.0, where the darker shading implies 
a higher number of events. 

to analyse their decays will be important to judge on the feasibility of the distinction based 
on the branching ratios. We postpone this study for future work. 



C. Differential Distributions in the Three-Body Decays 

Finally, we consider the differential distributions of the visible decay products in the 

three-body decays f^T + j + a/Gm terms of the quantity 

r(TR ^ T7i;a;^"S4"*) dx^dcosO ' ' ' ^ ' 

which is independent of the two-body decay, the total NLSP decay rate, and the Peccei- 
Quinn/Planck scale. In Fig. 5, the normalized differential distributions (19) with x™* = 
x^"* = 0.1 are shown for the axino LSP scenario (left) and the gravitino LSP scenario (right) 
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for m~ = 100 GeV, ms = 110 GeV, m^/m^ < 1, and mp^ = lOMeV.^ In the case of the 
gravitino LSP, the events are peaked only in the region where the photons are soft and the 
photon and the tau are emitted with a small opening angle ~ 0). In contrast, in the axino 
LSP scenario, the events are also peaked in the region where the photon energy is large and 
the photon and the tau are emitted back-to-back ~ tt). Thus, if the observed number 
of events peaks in both regions, there is strong evidence for the axino LSP and against the 
gravitino LSP. 

To be specific, with 10^ analysed stau NLSP decays, we expect about 165±13 (stat.) 
events for the axino LSP and about lOOilO (stat.) events for the gravitino LSP, which will 
be distributed over the corresponding ( X ■y J COS ^^)-planes shown in Fig. 5. In particular, in 
the region of > 0.8 and cos^^ < —0.3, we expect about 28% of the 165±13 (stat.) events 
in the axino LSP case and about 1% of the 100±10 (stat.) events in the gravitino LSP case. 
These numbers illustrate that O(IO^) of analysed stau NLSP decays could be sufficient 
for the distinction based on the differential distributions. To estabhsh the feasibihty of 
this distinction, a dedicated study taking into account the details of the detectors and the 
additional massive material will be crucial, which we leave for future studies. 

Some comments are in order. The differences between the two scenarios shown in Figs. 4 
and 5 become smaller for larger values of m^/m~. This ratio, however, remains close to unity 
for the stau NLSP in unified models. Furthermore, if mg < m~ < mLosp — mentioned as 
case (iv) in the Introduction — and r(f — > aX) ^ T{t GX), one would still find the 
distribution shown in the left panel of Fig. 5. The axino would then eventually decay into the 
gravitino LSP and the axion. Conversely, the distribution shown in the right panel of Fig. 5 
would be obtained if m~ < < mLosp — mentioned as case (iii) in the Introduction — 
and r(f aX) <^ r(f GX). Then it would be the gravitino that would eventually 
decay into the axino LSP and the axion. Barring these caveats, the signatures shown in 
Figs. 4 and 5 will provide a clear distinction between the axino LSP and the gravitino LSP 
scenarios. 

^ A similar comparison between the gravitino and a hypothetical spin-1/2 fermion with extremely weak 

Yukawa couplings was performed in Refs. [25]. Note that oTir result for the axino shown in Fig. 5 differs 
also from the one for the hypothetical spin-1/2 fermion due to different couplings. 
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V. CONCLUSION 



Assuming that a charged slepton is the NLSP, we have discussed signatures of both the 
gravitino LSP scenario and the axino LSP scenario in the framework of hadronic, or KSVZ, 
axion models [26]. These signatures can be observed at future coUiders if the planned 
detectors are equipped with 1-10 kt of additional material to stop and collect charged 
NLSPs [23, 24]. With calorimetric and tracking performance, this additional material will 
serve simultaneously as a real-time detector, allowing an analysis of the decays of the trapped 
NLSPs with high efficiency [23]. 

In the scenario in which the axino is the LSP, we have shown that the NLSP lifetime 
can be used to estimate the Peccei-Quinn scale /„. Indeed, if the axino is the LSP, the 
NLSP decays provide us with a new way to probe the Peccei-Quinn sector. This method 
is complementary to the existing and planned axion search experiments. The decays of the 
NLSP into the axino LSP will also allow us to determine the axino mass kinematically if it is 
not much smaller than the mass of the NLSP. The determination of both the Peccei-Quinn 
scale fa and the axino mass m~ will be crucial for insights into the cosmological relevance 
of the axino LSP. Once /„ and m~ are known, we will be able to decide if axinos are present 
as cold dark matter in our Universe. 

In the gravitino LSP scenario, the measurement of the stau NLSP lifetime can be used 
to determine the gravitino mass once the mass of the NLSP is known. This will be 
crucial for insights into the SUSY breaking mechanism. Moreover, if the gravitino mass 
can be determined independently from the kinematics and if the NLSP mass is known, the 
NSLP hfetime provides a microscopic measurement of the Planck scale [25]. Indeed, if the 
gravitino is the LSP, the hfetime of the NLSP depends strongly on the Planck scale and the 
masses of the NLSP and the gravitino. 

We have addressed the question of how to distinguish between the axino LSP and the 
gravitino LSP scenarios at colliders. If the mass of the LSP cannot be measured and if 
the NLSP lifetime is within the range 0(0.01 sec)-O(10 h), we have found that the NLSP 
lifetime alone will not allow us to distinguish clearly between the axino LSP and the gravitino 
LSP scenarios. The situation is considerably improved when one considers the three-body 
decay of a charged slepton NLSP into the associated charged lepton, a photon, and the LSP. 
We have found qualitative and quantitative differences between the branching ratios of the 
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integrated three-body decay rate with cuts on the photon energy and the angle between the 
lepton and photon directions. In addition, the differential distributions of the decay products 
in the three-body decays provide characteristic fingerprints. For a clear distinction between 
the axino LSP and the gravitino LSP scenarios based on the three-body decay events, at 
least of (9(10'^) of analysed stau NLSP decays are needed. If the mass of the stau NLSP is 
not significantly larger than 100 GeV, this number could be obtained at both the LHC and 
the ILC with 1-10 kt of massive additional material around the main detectors. 
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